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Chemical  Spill  Characteristics 
in  the  San  Diego  Bay 


ABSTRACT 


45  Dispersion  of  ocean  pollutants  in  estuarine  environments  and  bays  (such  as  San 

46  Diego  Bay)  depends  on  the  location  of  the  source  of  the  pollutants  relative  to  the 

47  mouth  and  the  tidal  excursion,  which  is  the  net  horizontal  distance  over  which  a 

48  pollutant  particle  moves  during  one  tidal  cycle  of  flood  and  ebb.  Pollutant  dispersion 

49  was  investigated  using  a  coupled  hydrodynamic  and  chemical  discharge  model  in 
so  this  study.  The  results  show  the  existence  of  two  distinct  (northern  and  southern) 

51  spill  patterns  of  pollutant  dispersion.  The  northern  spill  pattern  is  characterized  by 

52  fast  reduction  of  the  pollutant  concentration  in  the  water  column,  rapid  dispersion  of 

53  pollutants  to  the  San  Diego  port  and  to  outside  of  the  San  Diego  Bay,  and  slow  dis- 

54  persion  of  pollutants  to  the  southern  bay.  The  southern  spill  pattern  is  characterized 

55  by  slow  reduction  of  the  pollutant  concentration  in  the  water  column,  slow  disper¬ 
se  sion,  and  confinement  of  pollutants  in  the  southern  San  Diego  Bay.  The  results  may 

57  be  useful  for  ocean  pollution  control  and  management. 

58  Keywords:  Two  chemical  spill  patterns,  San  Diego  Bay,  ocean  pollution,  water 

59  quality  management,  chemical  dispersion 
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1.  Introduction 

he  San  Diego  Bay,  located  at  the 
west  coast  of  southern  California, 
connects  to  the  Pacific  Ocean  through 
a  single  channel  at  the  mouth  (Fig¬ 
ure  1).  It  is  a  semienclosed  bay  and  a 
natural  harbor  sheltered  by  overlapping 
peninsulas  (in  the  west,  Point  Loma, 
and  in  the  east,  Coronado).  The  bay 
has  been  intensively  engineered  to  ac¬ 
commodate  shipping  activities.  Ninety 
percent  of  all  available  marsh  lands  and 
50%  of  all  available  intertidal  lands 
have  been  reclaimed,  and  dredging  ac¬ 
tivities  within  the  bay  have  been  equally 
extensive  (Peeling,  1975).  The  shore¬ 
line  of  San  Diego  Bay  is  spotted  with 
high  pollution  from  shipbuilding  and 
ship  repair  facilities.  The  body  of  water 
in  the  bay  is  particularly  at  risk  because 
of  the  military  and  industrial  activities  in 
and  around  it.  Investigation  of  the  dis¬ 
persion  of  floating  chemicals,  such  as 
benzene,  is  very  important  for  the  mon¬ 
itoring  and  control  of  water  quality. 

The  San  Diego  Bay  has  a  “flipped  T” 
shape  and  is  nearly  25  km  long  and 
1-4  km  wide  (Figure  la).  The  bot¬ 
tom  topography  of  the  bay  is  not  ho¬ 
mogeneous,  with  an  average  depth  of 
6.5  m.  The  northern/outer  part  of 
the  bay  is  narrower  (1—2  km  wide) 
and  deeper  (reaching  a  depth  of 
15m),  and  the  southern/inner  part  is 


60  wider  (2—4  km  wide)  and  shallower 

61  (depth  less  than  5  m)  (Figure  lb). 

62  Once  pollutants  are  released  into  the 

63  San  Diego  Bay,  dispersion  of  pollu- 
64tants  depends  upon  the  hydrody- 
65namic  forcing  caused  by  exchange 
66  between  the  San  Diego  Bay  and  the 
67Pacific  Ocean  through  a  single 

68  north-south  channel,  which  is  about 

69  1.2  km  wide,  bounded  by  Point 

70  Loma  to  the  west  and  Zuniga  jetty 

71  to  the  east,  with  depths  between  5 

72  and  15  m.  The  west  side  of  the  chan- 

73  nel  is  shallower  than  the  east  side. 

74  Such  topographic  features  cause  a 
75 phenomenon  called  “tidal  pumping,” 

76  due  to  the  asymmetry  between  the 

77  flow  during  the  ebb  and  flood  tides 
78(Fischer  et  ah,  1979).  Transport 
79  time  for  pollutant  particles  moving 
so  out  of  the  bay  depends  on  the  hori- 
sizontal  distance  relative  to  the  mouth 

82  and  the  tidal  excursion,  which  is  the 

83  net  horizontal  distance  over  which  a 


water  particle  moves  during  one 
tidal  cycle  of  flood  and  ebb.  Numeri¬ 
cal  modeling  and  chemical/isotopic 
tracer  analyses  are  generally  used  to 
investigate  such  dependence  for 
water  quality  control  and  manage¬ 
ment.  Between  them,  numerical 
modeling  is  cost-effective  without  af¬ 
fecting  the  water  environment.  Here, 
a  numerical  modeling  study  is  pre¬ 
sented.  The  model  has  hydrodynamic 
and  chemical  discharge  components. 
The  hydrodynamic  part  is  driven  by 
tides  and  winds  and  predicts  the  ve¬ 
locity  field.  The  chemical  discharge 
part  is  driven  by  the  velocity  field 
from  the  hydrodynamic  model  and 
predicts  the  pollutant  dispersion. 

2.  Background 

2.1.  Vertically  Well-Mixed  Basin 

The  Space  and  Naval  Warfare  Sys¬ 
tems  Command  (SPAWAR)  deployed 
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FIGURE  1 


San  Diego  Bay:  (a)  main  geographical  locations  and  (b)  bathymetry. 
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three  acoustic  Doppler  current  profil¬ 
ers  (ADCPs)  in  the  San  Diego  Bay 
in  1993  (Figure  2)  with  a  broadband 
ADCP  (station  bb)  located  at  the 
mouth  of  the  bay  (32°42/25.8"N, 
1 17°13r30.6"W)  from  June  22  to 
July  23,  and  two  narrowband  ADCPs 
inside  the  bay:  station  nbl  located  at 
(32°'43.98"N,  1 17°12'55.68"W) 
from  June  22  to  August  26  and  sta¬ 
tion  nb2  located  at  (32°42r17.22"N, 
1 17°10/8.88"W)  from  June  23  to 


ns  August  27.  Figure  3  shows  time  series 

119  of  horizontal  velocity  components 

120  ( u ,  v)  at  three  different  depths  (surface, 

121  middepth,  and  bottom)  of  two  ADCP 

122  stations  (nb  1  and  nb2)  inside  the  bay. 

123  The  three  curves  are  very  close  together 

124  for  each  component  (u  or  v)  at  each 
125 station  (nbl  or  nb2),  showing  well- 
remixed  characteristics.  The  correlation 
127  coefficient  between  the  surface  and 
i28bottom  currents  is  97.2%  for  the 
129 u  component  and  96.3%  for  the 


v  component  at  station  nbl  and  92.0% 
for  the  u  component  and  94.7%  for 
the  v  component  at  station  nb2. 

2.2.  Atmospheric  Conditions 

From  National  Oceanic  and  Atmo¬ 
spheric  Administration’s  (NOAA) 
weather  description,  wind  forcing  is  al¬ 
ways  less  significant  than  tidal  forcing 
in  the  San  Diego  Bay.  The  mean  west¬ 
erly  winds  in  the  afternoon  and  mean 
easterly  winds  in  the  evening  and 
morning  are  less  than  5  m/s  with  prac¬ 
tically  no  storms  in  June,  July,  and  Au¬ 
gust.  Rain  occurs  mostly  in  winter  and 
almost  never  in  summer,  with  an  an¬ 
nual  precipitation  of  about  0.26  m. 
In  terms  of  estuarine  classification, 
the  San  Diego  Bay  is  generally  positive, 
i.e.,  drainage  inflow  exceeds  evapora¬ 
tion  (Pritchard,  1952).  However,  dur¬ 
ing  the  summer,  the  evaporation  rate 
(about  0.16  m)  exceeds  precipitation 
(near  zero)  (Peeling,  1975),  and  a  “re¬ 
versed  estuary”  phenomenon  is  ob¬ 
served  (Defant,  1961).  Small  water 
mass  flux  at  the  surface  (mostly  in  win¬ 
ter)  and  weak  wind  forcing  make  the 
San  Diego  Bay  a  tidally  driven  basin 
(Fagherazzi  et  al.,  2003). 

2.3.  Water  Quality 

Military  and  civilian  vessel  activi¬ 
ties  provide  sources  of  the  toxicity. 
Widespread  toxicity  in  the  San  Diego 
Bay  sediments  contains  copper,  zinc, 
mercury,  polycyclic  aromatic  hydro¬ 
carbons,  polychlorinated  biphenyls, 
and  chlordane.  No  single  chemical  or 
chemical  group  has  a  dominant  role 
in  contributing  to  the  identified  toxic¬ 
ity.  The  semienclosed  Shelter  Island 
Yacht  Basin  (a  boat  harbor)  has  been 
added  to  California’s  list  of  impaired 
water  bodies.  The  toxicity  comes 
from  specially  formulated  paints  that 
are  impregnated  with  biocides  and 
applied  to  boat  hulls  to  retard  the 
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FIGURE  2 

Location  of  the  ADCP  stations  deployed  by  SPAWAR  in  June  to  August  1993.  Note  that  station 
bb  is  located  at  the  mouth  of  the  San  Diego  Bay. 


176  growth  of  fouling  organisms  such  as  isiU.S.  naval  bases,  the  San  Diego  Bay 

177  barnacles.  182  is  a  possible  target  of  chemical  attack 

178  In  the  current  environment  of  i83with  many  possible  chemical  com- 

179  threats  to  homeland  security  and  as  a  184  pounds.  For  example,  benzene  is  an 
iso  big  city  waterway  that  hosts  large  185  organic  chemical  compound  with  the 


FIGURE  3 


Time  series  of  (u,v)  components  from  station  nbl  at  surface  (yellow),  middle  depth  (purple), 
and  bottom  (blue)  for  station  nbl  (top)  and  nb2  (bottom):  (a)  [/component  and  (b)  v compo¬ 
nent.  (Color  versions  of  figures  available  online  at:  http://www.ingentaconnect.com/content/mts/ 
mtsj/201 1/00000045/00000002.) 
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molecular  formula  QH^.  It  is  some-  186 
times  abbreviated  Ph-H.  Benzene  is  a  187 
colorless  and  highly  flammable  liquid  188 
with  a  sweet  smell,  an  aromatic  hydro-  189 
carbon  and  the  second  [/fi-annulene  190 
([6]-annulene),  and  a  cyclic  hydrocar-  191 
bon  with  a  continuous  pi  bond.  It  is  192 
also  related  to  the  functional  group  193 
arene,  which  is  a  generalized  structure  194 
of  benzene.  Here,  we  use  benzene  as  an  195 
example  to  show  the  effect  of  tidal  196 
pumping  on  the  chemical  spill  patterns  197 
in  the  San  Diego  Bay.  Sewage  runoff  is  198 
important  but  not  included  in  this  199 
study.  200 


3.  Hydrodynamic  Chemical  201 
Discharge  Model  202 

3.1.  Water  Quality  Management  203 
and  Analysis  Package  204 

Water  Quality  Management  and  205 
Analysis  Package  (WQMAP)  is  a  nu-  206 
merical  hydrodynamic  model  devel-  207 
oped  at  Applied  Science  Associates,  208 
Inc.  (ASA)  with  fitted  boundaries  209 
(Muin  and  Spaulding,  1996,  1997).  210 

The  model  is  configured  to  run  in  a  211 
vertically  averaged  (barotropic)  mode  212 
or  as  a  fully  three-dimensional  (baro-  213 
clinic)  mode.  Several  assumptions  are  214 
made  in  the  model  formulation,  215 
including  hydrostatic  approxima-  216 
tion,  Boussinesq  approximation,  and  217 
incompressibility.  In  this  study,  the  218 
two-dimensional  version  is  used.  219 

WQMAP  was  implemented  for  the  220 
San  Diego  Bay,  covering  an  area  of  221 
43  km2.  The  computational  mesh  222 
has  150  x  200  (30,000)  grid  nodes  223 
with  an  average  horizontal  resolution  224 
of  40  m.  The  sources  for  the  water  225 
depths  are  the  NOAA  sounding  data  226 
and  navigation  charts  and  the  navy-  227 
conducted  bathymetry  survey.  The  228 
navy  data  shows  that  the  water  depths  229 
in  regions  near  the  bay  entrance  are  230 
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significantly  deeper  than  the  water 
depths  shown  on  the  NOAA  naviga¬ 
tion  chart  (Wang  et  ah,  1998).  The 
most  up-to-date  bathymetry  data  are 
used  in  the  model. 

The  model  was  span  up  from  a 
quiescent  initial  condition  and  uni¬ 
form  temperature  (16°C)  and  salinity 
(34  ppt)  for  1  day  and  then  integrated 
with  tidal  and  wind  forcing  from  time 
00:00  on  22  June  1993  to  23:54  on  27 
August  1993  with  time  step  of  6  min. 
The  CFL  condition  is  satisfied  at  this 
time  step.  Sea  surface  elevation  at  the 
mouth  of  San  Diego  Bay  is  available 
every  6  min  at  the  NOAA  Station 
9410170,  located  at  (32°42'48"N, 
117°10/24"W)  and  taken  as  the  tidal 
forcing  function.  The  integration 
period  is  selected  from  22  June  1993 
to  27  August  1993  (see  Figure  3)  in 
accordance  with  the  observational 
period  of  three  ADCPs  for  model- 
data  intercomparison. 


255  Statistical  analysis  (Chu  et  al. ,  200 1 ) 

256  shows  good  correlation  between  mod- 

257  eled  and  observed  horizontal  velocity 

258  with  the  correlation  coefficients  above 

259  0.90  in  all  cases.  At  nb  1 ,  the  correlation 

260  coefficient  of  the  u  component  is  0.92. 

261  The  observational  u  component  ranges 

262  between  -51.8  and  44.5  cm/s,  and  the 

263  modeled  u  component  changes  be- 

264  tween  -46.9  and  40.8  cm/s  (Figure  4). 
265The  correlation  coefficient  of  the 
266  v  component  is  also  0.92.  The  observa- 
267tional  v  component  ranges  between 

268  -3 1 .6  and  29.6  cm/s,  and  the  modeled 

269  v  component  changes  between  -37.0 

270  and  32.0  cm/s.  Overall,  the  model  ve- 

271  locities  are  reasonably  good,  especially 

272  taking  into  account  that  the  data  and 

273  the  model  output  are  not  at  exactly 

274  the  same  geographic  location  and  the 

275  proximity  of  the  ADCPs  to  shore.  If 

276  finer  grid  and  more  accurate  bathyme- 

277  try  are  used,  the  model  results  may  be 

278  further  improved. 


FIGURE  4 

Model  (blue  curve)  and  (ADCP)  data  (purple  curve)  comparison  for  station  nbl  (top )  and  nb2  (bottom): 
(a)  [/component  and  (b)  1/ component.  (Color  versions  of  figures  available  online  at:  (Color  ver¬ 
sions  of  figures  available  online  at:  http://www.ingentaconnect.com/content/mts/mtsj/2011/ 
00000045/00000002.) 


3.2.  Chemical  Discharge  Model  279 

A  chemical  discharge  model  (called  280 
CHEMMAP)  was  also  developed  at  28i 
ASA  to  predict  or  to  simulate  surface  282 
and  subsurface  spills,  slick  spreading,  283 
transport  of  floating,  dissolved  and  284 
particulate  materials,  evaporation  and  285 
volatilization,  dissolution  and  adsorp-  286 
tion,  sedimentation,  and  degradation.  287 
The  model  inputs  are  density,  viscos-  288 
ity,  vapor  pressure,  surface  tension,  289 
water  solubility,  environmental  degra-  290 
dation  rates,  and  adsorbed/dissolved  291 
partitioning  coefficients.  The  model  292 
outputs  are  the  trajectory  and  fate  293 
of  floating,  sinking,  evaporating,  294 
soluble/insoluble  chemicals,  and  esti-  295 
mation  of  the  distribution  of  chemical  296 
elements  (mass  or  concentration)  on  297 
the  surface,  in  the  water  column,  and  298 
in  the  sediments.  The  model  separately  299 
tracks  surface  slicks,  entrained  droplets  300 
or  particles  of  pure  chemical,  chemical  301 
adsorbed  to  suspended  particulates,  302 
and  dissolved  chemicals  (McCay  and  303 
Isaji,  2002).  More  specifically,  the  304 <33 
model  can  predict  the  swept  area  by  a  305 
floating  chemical,  as  well  as  total,  ab-  306 
sorbed,  dissolved,  and  particulate  con-  307 
centration  in  both  the  water  column  and  308 
sediments,  and  can  determine  the  range  309 
and  direction  of  contamination  caused  310 
by  the  spill  at  a  particular  location.  311 

4.  Chemical  Spill  Patterns  312 

Suppose  that  one  barrel  of  a  313 
chemical  (e.g.,  10  tons  of  benzene)  is  314 
released  into  the  water  from  a  small  315 
boat  at  00:00  on  day  1  at  (1)  north-  316 
ern  San  Diego  Bay  (32°43'N,  117°  317 

13.05'W)  (point  2  in  Figure  la)  and  318 
(2)  southern  San  Diego  Bay  (32°39/N,  319 

117°07.92'W)  (point  4  in  Figure  la).  320 
The  release  depth  is  1  m,  and  the  initial  321 
plum  thickness  is  0.5  m.  Two  distinct  322 
spill  (northern  and  southern)  patterns  323 
are  found  for  all  the  chemicals.  Here,  324 
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spill  patterns  of  benzene  are  presented 
for  illustration. 


FIGURE  5 

Benzene  dissolved  concentration  out  of  the  bay  1 2  h  after  being  dropped  in  the  North  San  Diego  Bay. 


4.1.  Northern  Spill  Pattern 

After  the  pollutants  are  released  at 
the  northern  San  Diego  Bay  (32°43/N, 
1 17°13.05'W),  the  pollutants  disperse 
generally  from  the  northern  bay  (north 
of  32°39'N)  to  outside  of  the  San 
Diego  Bay.  They  disperse  very  little 
into  the  southern  bay  (south  of  32° 
39  rN).  The  benzene  reaches  the  San 
Diego  port  (Figure  la)  in  about  3  h. 
It  transports  outside  of  the  San  Diego 
Bay  in  12  h  (Figure  5).  The  southern 
bay  is  not  contaminated  for  the  first 
5  days  (Figure  6a)  and  weakly  affected 
after  32  days  (Figure  6b).  Rapidly 
weakening  of  the  pollutant  concentra¬ 
tion  in  the  water  column  is  found.  The 
pollutant  concentration  is  20%  after 
5  days,  reduces  to  10%  after  15  days, 
and  reaches  4%  after  30  days  (Fig¬ 
ure  7).  There  is  plenty  of  time  to  take 
protective  measures  for  the  southern 
bay  (Chula  Vista  area),  where  the  impact 
of  such  an  incident  would  be  minor. 

4.2.  Southern  Pattern 

After  the  pollutants  are  released  at 
southern  San  Diego  Bay  (32°39/N, 
117°07.92'W),  the  spill  pattern  is  to¬ 
tally  different  from  the  northern  spill 
pattern.  The  pollutants  disperse  gen¬ 
erally  inside  the  bay  with  very  few  pol¬ 
lutants  reaching  the  32°4lrN  parallel. 
However,  the  naval  station  (Figure  la) 
is  affected  within  12  h  (Figure  8a)  and 
completely  contaminated  in  less  than 
3  days.  It  is  important  for  protective 
measures  to  highlight  this  pattern 
because  a  chemical  attack  in  the  south¬ 
ern  part  of  the  bay  would  affect  the 
naval  station.  After  17  days,  the  dis¬ 
solved  benzene  reaches  the  San  Diego 
port  (Figure  8b).  After  32  days,  the  dis¬ 
solved  benzene  is  confined  in  the 
southern  San  Diego  Bay  (Figure  9).  It 


FIGURE  6 

Dispersion  of  benzene  (a)  5  days  and  (b)  32  days  after  being  dropped  in  the  North  San  Diego  Bay. 
(A)  (B) 


FIGURE  7 

Mass  balance  for  benzene  dropped  in  the  northern  San  Diego  Bay. 
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FIGURE  8 


379 


Benzene  concentration  (a)  12  h  and  (b)  17  days  after  being  dropped  in  the  South  San  Diego  Bay. 


FIGURE  9 

Dispersion  of  benzene  32  days  after  being  dropped  in  the  southern  San  Diego  Bay. 


clearly  shows  that  the  pollutants  are 
more  likely  confined  in  the  southern  380 
San  Diego  Bay  for  quite  a  long  period.  38i 
Temporal  variability  of  the  pollutant  382 
concentration  in  the  water  column  is  383 
quite  different  between  the  southern  384 
(Figure  10)  and  northern  (Figure  7)  385 

spill  patterns.  Slow  reduction  of  the  386 
pollutant  concentration  in  the  water  387 
column  is  found  for  the  southern  388 
spill  pattern.  The  pollutant  concentra-  389 
tion  is  more  than  30%  after  10  days,  390 
reduces  to  25%  after  15  days,  and  391 
reaches  10%  after  30  days  (Figure  10).  392 

This  pattern  may  affect  human  beings  393 
and  the  environment  as  a  result  of  the  394 
longer  period  of  confinement  of  pollu-  395 
tants  in  the  southern  San  Diego  Bay.  396 

5.  Conclusions  397 

In  this  study,  two  distinct  (north-  398 
ern  and  southern)  chemical  spill  pat-  399 
terns  were  found  depending  on  the  400 
location  of  the  pollutant  source.  The  401 
northern  spill  pattern  occurs  when  402 
the  pollutants  are  released  in  the  north-  403 
ern  San  Diego  Bay.  It  is  characterized  404 
by  fast  reduction  of  the  pollutant  con-  405 
centration  in  the  water  column,  rapid  406 
dispersion  of  pollutants  to  the  San  407 
Diego  port  and  to  outside  of  the  San  408 
Diego  Bay,  and  slow  dispersion  of  pol-  409 
lutants  to  the  southern  bay.  The  south-  410 
ern  spill  pattern  appears  when  the  411 
pollutants  are  released  in  the  southern  412 
San  Diego  Bay.  The  southern  spill  pat-  413 
tern  is  characterized  by  slow  reduction  414 
of  the  pollutant  concentration  in  the  415 
water  column,  slow  dispersion,  and  416 
confinement  of  pollutants  in  the  417 
southern  San  Diego  Bay.  Although  418 
the  modeling  results  are  useful,  one  419 
should  be  precautious  in  applying  420 
them  to  ocean  pollution  monitoring,  421 
control,  and  management.  This  is  due  422 
to  uncertainties  in  the  numerical  model  423 
such  as  the  bathymetry,  discretization,  424 
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FIGURE  10 

Mass  balance  for  benzene  dropped  in  the  southern  San  Diego  Bay. 
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boundary  configuration,  and  forcing 
functions.  Another  problem  is  the 
lack  of  recent  data  for  the  San  Diego 
Bay.  The  comparison  was  conducted 
between  hydrodynamic  model  output 
and  old  ADCP  observations  because  of 
the  lack  of  more  recent  data.  These  is¬ 
sues  need  to  be  carefully  considered  be¬ 
fore  using  these  results. 
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